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ABSTRACT 

Two popular models of optically thick relativistic outflows exist: the wind and the shell. We propose a unified 
treatment of photospheric emission within these models. We show that quite counterintuitive situations may 
appear when e.g. geometrically thin shell may behave as thick wind. For this reason we introduce notions of 
photon thick and photon thin outflows. They appear more general and better physically motivated than winds 
and shells when photospheric emission is considered. 

We obtain light curves and observed spectra for both photon thick and photon thin outflows. In the photon 
thick case we generalize the results obtained for steady wind. It is our main finding that the photospheric 
emission from the photon thin outflow is dominated by diffusion and produces non thermal time integrated 
spectra, which may be described by the Band function well known in the GRB literature. 

Energetic GRBs should produce photon thin outflows and therefore when only time integrated spectra for 
such GRBs are available we naturally expect them to have a Band shape. In the literature Band spectra for the 
photospheric emission of GRBs are obtained only involving additional dissipative mechanisms which are not 
required here. 

Subject headings: Opacity — Plasmas — Radiation mechanisms: thermal — X-rays: bursts 



1. INTRODUCTION 

High Lorentz factors of the bulk motion of various outflows 
are common in relativistic astrophysics. The best known ex- 
amples are Active Galactic Nuclei (Maraschi 2003), micro- 
quasars and Gamma-Ray Bursts (GRBs) (.Piraa.2004.) . In the 
latter case outflows indeed reach ultrarelativistic velocities. 

Various models are suggested to explain the acceleration of 
outflows to ultrarelativistic velocities. The electromagnetic 
model (Lyutikov 2006) assumes that the energy in the source 
of GRB is converted into electromagnetic energy which is 
transported in the form of a Poynting flux. We adopt here an- 
other popular idea that the energy release leads to creation of 
an optically thick source which expands due to thermal accel- 
eration. This idea is the basis of both the fireball (Piran ( 1999", 
[2004) and references therein) and the fireshell (Ruffini et al., 
yOOl, 2009) and refer ences therein) mod els. 

In pioneer works by 'Goodman' ("1986^) who considered an 
instant explosion, and by Paczynski (1986) who discussed 
a gradual energy release, a conclusion was reached that the 
electron-positron plasma is created in the source of GRB. As- 
suming further that the plasma reaches thermal equilibrium 
they focused on hydrodynamic expansion in such models and 
gave photometric and spectroscopic predictions for GRBs. 
Later, baryonic loading of fireballs was considered for explo- 
sions by Shemi & Piran (1990) and for winds by Paczynski 
(|J990). Abramowicz et al.l (119 91) considered the appearance 
of the photosphere of the relativistic wind to a distant observer 
showing that its shape is concave. 

The interest in photospheric emission from rela- 
tivistic winds has bee n revive d re cently in pap ers b y 
Daiane & Mochkovitch" ('2002'), 'Pe'eretal.' ('2007), 
Beloborodov (2011), Ryde et al. (2011a), Pe'er&Ryde 
(1201 ID . and others. In the fireshell model, which assumes 
an explosive energy release, the first potentially visible 
component of any GRB, the Proper GRB (iBianco et al.ll2001t 
iRuffini et al.ll200 T). comes from the thermal flash of radiation 
emitted when the outflow becomes transparent for photons. 



One of the main scopes of this paper is to show that the as- 
soci ation of "instantan eous energy release" with "thin shell" 
(e.g. iPiran et all (1 1 9931) and Meszar oTetaH (119931)) and "con- 
tinuous energy release" with "thick wind" (e.g. iPaczvnskil 
( 1986, 1990)) generally adopted in the literature is incomplete 
with respect to the photospheric emission. A unified treat- 
ment of ultrarelativistic outflows originating both from instan- 
taneous and from continuous energy release is presented here. 
We also propose a new classification of such outflows which is 
complementary to the traditional division on shells and winds. 
When photospheric emission is considered this new classi- 
fication provides some additional physical insights. Based 
on these results we present a semi-analytic treatment for ob- 
served flux of photospheric emission assuming isotropic ther- 
mal distribution of scattered photons in the comoving frame. 
We compute both instantaneous and time-integrated spectra. 
These results are then applied within both the shell and the 
wind models of GRBs. Remarkably the time-integrated spec- 
trum of energetic GRBs is predicted to have a Band shape 
(iBand et al.lll993h . 

The structure of the paper is as follows. In Section 2 we dis- 
cuss, compare and contrast an impulsive explosion and grad- 
ual energy release, giving rise respectively to an ultrarelativis- 
tic shell and wind. In Section 3 we present the computation 
of the optical depth and photospheric radius of the relativistic 
outflow and discuss applicability of recovered asymptotics to 
GRBs. In Section 4 we discuss radiation diffusion in an ex- 
panding relativistic outflow. Observed light curves and spec- 
tra of photon thick and photon thin outflows are computed in 
Sections 5 and 6, respectively. In Section 7 we discuss main 
results obtained in this work and their implications for GRBs. 
Conclusions follow. Details of derivation and application of 
the radiative diffusion approximation are given in Appendix. 



2. OPTICALLY THICK RELATIVISTIC OUTFLOWS: 
WIND VS. EXPLOSION 
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Relativistic outflows are generafly classified as winds or 
shells depending on how fast the energy in their sources is 
released. Consider energy release in a source of linear dimen- 
sion /?(). If the time scale of energy release is short At ~ Rq / c, 
an explosion occurs, which may be characterized by the size 
Ro, total energy released £0 and total baryonic mass M. Rela- 
tivistically expanding material forms a shell having width ap- 
proximately When the energy is released gradually, on a 
time scale At ^ Rq /c, but the source luminosity L exceeds 
the Eddington limit, a wind is formed, which is characterized 
by its activity time Af , luminosity L and mass ejection rate M. 

In both wind and shell models of cosmological GRBs the 
region of energy release Rq is clearly macroscopic, being of 
the order of the size of compact astrophysical objects, about 
Ro ~ 10** cm. Considering the isotropic energy involved £0 ^ 
lO^'* erg the temperature of plasma in this region is typically in 
MeV range. The time scale of thermalization in such plasma 
is of the order of lO'^^ sec (lAksenov et al.ll2007l I201Q) , be- 
ing much shorter than the dynamical time scale ^0 /c- Since 
such dense plasma is optically thick thermal equilibrium is 
established prior to expansion in both shell and wind models. 
Opacity in such plasma is dominated by Compton scattering. 
When the temperature in the source of the relativistic outflow 
is large enough for electron-positron pair creation, e'^e~ pairs 
make an additional contribution to the opacity. 

Whether the outflow becomes relativistic or not depends 
on the entropy in the region where the energy is released. 
Both the wind and explosion cases can be parameterized 
dShemi & PiranI [1990 ) by a dimensionless entropy parame- 
ter 7], equal to Eq/Mc^ for shell model and L/Mc^ for wind 
model. When the baryonic loading B = rf^ (iRuffini et al.l 
I2OOOI) is sufficiently small the baryons will be accelerated to a 
relativistic velocity v of bulk motion, attaining large Lorentz 



factors up to r = 



i-{v/cy 



-1/2 



T], while in the opposite 



case of large baryonic loading the outflow remains nonrela- 
tivistic with v ~ c/ \/2rj. 

In what follows we consider only ultrarelativistic spheri- 
cally symmetric outflows with F ^ 1 . In the simplest cases of 
wind or explosion in vacuum, the dynamics of the out flow is 
divide d into an acceleration phase and a coasting phase dPiranl 
Il999h with respectively 

r fry' 

Ro' """""V^j ' ^o<r<T]Ro: (1) 

r~ 77 = const, nc~ — [^] , r>i]Ro, (2) 

where tic is the comoving number density of baryons in the 
outflo'wQ. Notice that in the case of an impulsive explosion 
for r ^ Ro the matter and energy appear to a distant observer 
to be concentrated in a geometrically thi n shell having wi dth 
I ^ Ro due to the relativistic contraction (iPiran et al.llI993h . 

It is important to stress that both an infinitely long wind 
with a time-independent mass ejection rate and luminosity on 
the one hand, and an infinitely thin shell originating from an 
instantaneous explosion in infinitely thin region represent two 
limiting cases for the energy release. 

During both acceleration and coasting phases the continuity 

' All quantities with subscript "c" are measured in comoving reference 
frame, and all quantities without this subscript are measured in laboratory 
reference frame. 



equation for the laboratory number density reduces to n cx r ^. 
We take for the laboratory density profile 



no 



y I , R{t)<r <R(t) + l, 



(3) 



0, 



otherwise, 



where R{t) is the radial position of the inner boundary of the 
outflow. Such an outflow may be produced by a gradual en- 
ergy release with constant luminosity and mass ejection rate 
on a finite time Af and we will refer to it as the portion of 
wind. 

Below we show that both the shell and the wind defined 
above may appear for photons emitted inside it as long wind 
or as thin shell, depending on the initial conditions that spec- 
ify respectively their width / and activity duration Af. This 
appearance is a consequence of ultrarelativistic expansion of 
the outflow. It is crucial to keep in mind that photons emitted 
inside the expanding outflow propagate in a medium whose 
laboratory number density depends both on radial coordinate 
and on time n{r,t). For photons propagating in the wind the 
spatial dependence of the number density plays the key role, 
while for photons propagating in the shell its time dependence 
is crucial. 

3. OPTICAL DEPTH AND PHOTOSPHERIC RADIUS 
The optical depth along the photon world line £ is defined 



as 



(4) 



where a is Thomson cross section, /'^ is the 4-current of par- 
ticles, on which the photon scatters, and c/x^ is the element of 
the photon world line. 

Consider a spherically symmetric expanding outflow with 
an ultrarelativistic velocity v = /3c ~ 1-1 /2r^. Assume that 
the photon is emitted at time t from the interior boundary r = R 
of the outflow and it pr opagates outwards. This assumption is 
relaxed in Section ISTI The optical depth computed along the 
photon path from (|4]i is (see e.g. lAbramowicz et al.l (1199 II) ) 



R+AR 



CT« (1-/3 COS 61) 



dr 



cost 



(5) 



where R + AR is the radial coordinate at which the photon 
leaves the outflow, and 6 is the angle between the velocity 
vector of the outflow and the direction of propagation of the 
photon, n is the laboratory number density of electrons and 
positrons, which may be present due to pair production. Ge- 
ometry of the outflow and used variables are illustrated by Fig. 
□ 

A pure electron-positron pl asma rea ches thermal equilib- 
rium before expansion ( Aksen ov et al.l l2007. 2010) and it re- 
mains accelerating until it becomes transparent to radiation. 
Due to exponential dependence of thermal pair density on the 
radial coordinate transparency is reached at 



kT± ~ 0.040m,c^ 



(6) 



rather independent of the initial conditions. Note that the op- 
tical depth f or an expanding electr on-positron-photon shell 
computed by IShemi & PiranI ( [T990I) is incorrect since it uses 
photon thin asymptotics, see 3rd line in d?) below, which 
never applies to the pure e^e~ outflows. The formula I© is 
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Fig. 1 . — Geometry of the outflow and variables used. Observer is located 
to the right at infinity. 



in agreement with works of iGrimsrud & WassermanI (Il998h 
and|Lr&Sm(2008). 

Electron-positron plasma with baryonic loading also 
reaches thermal equilibr ium before its expansion starts 
jAksenov et al.l2009ll2010 ). With decreasing entropy rj opac- 
ity due to electrons associated with baryons increases and 
eventually dominates over pair opacity. Straightforward cal- 
culations with the laboratory density profile gives 



R()l o 



(7) 



where 



To = anoRo = 



aL 



AirmpC^Rolr] Am^ipC^Rorj 



(8) 



Let us interpret the formula dTji. On the one hand, first two 
lines in (jTjl correspond to the case when the photon propagates 
inside the outflow for a significant time so that the number 
density on its path substantially decreases before it leaves. In 
this respect the outflow is "a long wind", even if the labora- 
tory thickness of the outflow may be small, / <C R- We refer 
to this case as a photon thick outflow. On the other hand, the 
third line in d?) corresponds to the case when the number den- 
sity of the outflow does not change substantially on the photon 
world line before it escapes. In this respect the outflow is "a 
thin shell" even if the duration of explosion could be long and 
a wind was launched. We refer to this latter case as a photon 
thin outflow. For instance, a geometrically thin ultrarelativis- 
tically expanding shell may be both thin or thick with respect 
to the photon propagating inside it. 

Similar consideration may be applied to a photon emitted 
at any distance ^ from the outer boundary of the outflow, see 
Section ISTI It is clear then, that even in a photon thick out- 
flow there is always a photon thin layer located near the outer 
boundary. During acceleration phase such a photon thin part 
accounts for a fraction not larger than F"' of the entire width 
of the outflow^. 

In the derivation of ^ we considered a portion of a rela- 
tivistic wind, but these results are generic and apply to any 
density profile of the outflow. 

^ Formally photon thin accelerating solution exists, and it is given by the 
last line in Q. However, its validity condition is / Sq/R = So /r. 



The photospheric radius Rpi, is defined by equating (|7]i to 
unity. It is worth noting that it is the lower limit in the integral 
(O that is associated with the photospheric radius, but not the 
upper one. The upper limit in Q is the radius at which the 
photon leaves the outflow, even if it may decouple from the 
outflow at much smaller radius, as in the photon thick case. 
The characteristic parameter which determines the type of the 
outflow at the photospheric radius is the ratio tq /(4F"^). When 
it is much smaller than one the outflow is photon thick, while 
in the opposite case the outflow is photon thin. 

Consider now the shell model of GRBs with typical pa- 
rameters expressing their total energy as £0 = 10^'*£'54 erg, 
initial size as Rq = I = 10^ Ra cm and entropy parameter as 
7] = We find the following asymptotic solutions for the 

photospheric radius together with domains of their applicabil- 
ity 



4.4 X 10"'(£54«8)'''^ cm, 



£54<4.8x 10- 



-20, 



Rph = S 



1.8 X 10'2(£547y-i/;8)'''^ cm, 

4.8 X 10-^%^R^^ < £54 < 3.2 X 10-^4rI 

1.8 X 10"E54TJ2^R^^ cm, 
3.2 X 10-^r]^Rl < £54 < 1-1 X lO-^v^Rl 

1/2 



(9) 



5.9 X 10"* (£5477-1)''" cm, 



£54 > 1.1 X lO-^j]^Rf 



For very small baryonic loading, or in other words for a 
pure electron-positron plasma, the photospheric radius does 
not depend on 77 parameter For increasing baryonic loading it 
increases as 77"'''"' (accelerating photon thick solution). In both 
these cases the Lorentz factor at the photosphere is not equal 
to 1], but it is much smaller. For larger baryonic loading the 
photospheric radius steeply increases as 77"^ (coasting photon 
thick solution), and finally it increases as 77"'/^ (coasting pho- 
ton thin solution), see Fig. |2] There we also show as function 
of the entropy parameter the following quantities computed at 
the photospheric radius: the Lorentz factor, the observed and 
comoving temperatures, and fraction of energy emitted from 
the photosphere to the total energy, for different values of the 
total energy Eq. Dependence of parameters of transparency 
on initial conditions is also illustrated in Tab. [T] It is clear that 
the highest Lorentz factors at photospheric radius are attained 
in photon thick asymptotics. The largest transparency radii 
are reached instead in photon thin asymptotics. 

Fig. [3] shows the energy-baryonic loading diagram, where 
the regions of validity of the asymptotics discussed above are 
indicated explicitly for typical parameters of GRBs. For all 
the relevant range of GRBs parameters 10"*^ erg < £0 < 10^^ 
erg and 10^ cm < /?() < 10'^ cm all four asymptotics are 
present in the interval 10 < 77 < lO'". 
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TABLE 1 

Dependencies of transparency parameters on initial conditions: energy Eg, entropy tj, and radius Rq for the shell model. 



Regime of transparency 


Rpl, 


r,,;, 


kTpi, 




Pair 

Acceleration 
Coasting photon tliick 
Coasting photon thin 


£l/3 1/3 

Eo fio' 


£l/3 -1/3 j^-Vi 

V 


£-1/12 1/3 ^-1/12 

£-S/12 5/3 /jV/12 

-1/12 -1/12 
^0 '^O 


^1/4 ^-3/4 

£-5/12 8/3 ^7/12 
-1/12 -1/12 




Fig. 2. — From top to bottom: photospheric radius Rpi,, Lorentz factor F,,;,, 
observed Tots and comoving Tpi, temperatures, and ratio of the energy emitted 
to the total energy Epi,/Eo at photospheric radius as functions of entropy rj 
for shells with different total energy Eq but the same width 1 = Rq = IQ^ cm. 
All four regimes with different asymptotics are clearly visible and dashed 
black lines corresponding to their domain of validity from Eq. j9) are shown. 
Curves are drawn for £"0 equal to: 104' erg (green), 10^' erg (blue), 10^^ erg 
(violet), and 10^' erg (red). Dashed thick lines denote the diffusion radius for 
each energy. 
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Photon Thin 



£? 10 




Accelerating 
Photon Thick 



105^* 
10^2 
1050 



10' 



,48 



10^ 
V 

Fig. 3. — The energy-baryonic loading diagram showing the validity of the 
various asymptotic solutions for the photospheric radius for typical parame- 
ters of GRBs with / = = 10** cm. 

In the case of gradual energy release resulting in relativistic 
wind an additional parameter is present, that is the duration 
of energy release, which we parameterize as At = 1 Afi s. In- 
stead of the total energy £0 the luminosity L = lO^^Lso erg/s is 
used. The corresponding photospheric radius is 

' 8.1 X lO^^Lj^^Vg/^ cm, 



L5()<5.3x lO-'X^: 



1.3x 10i"(L5o?72'/?i)'^^ cm, 

5.3 X lO-^^4Rf < L50 < 9.8 X lO-^4Rsi, 

5.9 X 10i°L5or/2^ cm, 

9.8 X 10-^4Rs<^L5o < 10^/1 Afi, 



5.9x lO'OfLsoAfiTyj'V''^ cm, 



(10) 



L50 > lO^Ty^Afi 



In Fig. in we show as function of the entropy parameter the 
following quantities computed at the photospheric radius: the 
Lorentz factor, the observed and comoving temperatures, and 
fraction of energy emitted from the photosphere to the total 
energy, for different duration of the wind with the total en- 
ergy £■() =10^' erg, and inner boundary radius Rq = 10^ cm. 
Wind duration ranges from 10 ms to 10 s. The correspond- 
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TABLE 2 

Dependencies of transparency parameters on initial conditions: energy release rate L, entropy rj, duration Ar, and radius Rq 

FOR THE WIND MODEL. 



Regime of transparency 


Rpi, 








Pair 


Ll/4 




0, 






r,-l/3 






0.040mfc2 




Ll/4 


^-1/2 

, 
r,-l/2 
'^0 


Acceleration 












-1/12 


^1/3 




Ll/4 


Coasting plioton thick 


L 






V 




-5/12 


^5/3 


„l/6 
'^O 


L-V'2 ,,8/3 


r,l/6 
'^O 


Coasting plioton thin 


Ll/2 




Ari/2 




V 




L- 


-1/12 


Ari/3 




L-'/i2 ^ Ar'/' 


r,l/6 




Fig. 4. — The same as in Fig. |2]for winds with different duration, but 
the same total energy £"0 = 10^' erg and radius of origin So = 10** cm. All 
four regimes with different asymptotics are clearly visible. Curves are drawn 
for A; ranging from 10"~ s (green) to 10 s (red) in steps of one order of 
magnitude. 



ing wind luminosity varies from 10 erg/s to 10 erg/s. Fig. 
|5] shows the luminosity-baryonic loading diagram where the 
regions of validity of the asymptotics discussed above are in- 
dicated. Dependence of transparency parameters on the initial 
conditions of the wind is presented in Tab. |2] 
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Fig. 5. — The luminosity-baryonic loading diagram showing the validity of 
the various asymptotic solutions for photospheric radius of wind with dura- 
tion Ar = 0.1 s. Notation is the same as on Fig.ia 



3.1. Comparison with previous works 

The expressions for the opti cal depth of a r elativistic wind 
were obtained e.g. by Meszaros & ReesI (12000 ). Their formu- 
las coincide with our (|7]l up to a numerical factor which comes 
from the integration over the radial coordinate. It should be 
noted, however, that only the pho t on thi ck asymptotic limit 
is discussed in Meszaros & ReesI (12000 ). The photon thin 
asymptotic limit may also be valid for relativistic winds in 
the coasting phase, provided that / <^ R/{2rf-). This is an in- 
dependent condition from A? ^ Rq jc and it is therefore pos- 
sible to give the following constraints for Af under which the 
outflow takes the form of a wind, but it is photon thin at the 
photospheric radius: 



c 4r]^ c 



(11) 



photospheric radius for photon thick and photon thin asymp- 
totics for a portion of the coasting relat ivistic wind was ob- 
tained bv lDaigne & Mochkovitoiil (120021) . 

Similar considerations apply to an ultrarelativistic shell 
which is cons i dered e.g. in Shemi & Piran (1990) and 
'Mesz aros et"aD (1993') and by Ruffi ni et alJ (12000) in the pho- 
ton thin approximation. The corresponding condition that the 
shell at the photospheric radius appears to be photon thin is 
To ^ 477^. It is possible, however, that initial conditions sat- 
isfy the opposite constraint which results in a photon thick 
shell, see N akaretaP (120051) . 

Numerical hydrodynamic simulations produce complex 
density, temperature and Lorent z facto r pro files of the out- 
flow. In particular, iPiran et aTl (Il993h and iMeszaros et alj 
( 119931) considered an explosion in a compact region with ra- 
dius Rq and studied numerically the hydrodynamic evolution 
of an optically thick plasma with various initial conditions. 
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They have shown that the plasma forms a relativistically ex- 
panding shell with some density and velocity profiles. The 
characteristic width of the density profile appears to be con- 
stant up to large radii, but later it increases linearly with radius 
due to the fact that the Lorentz factor appears to be monoton- 
icaUy increasing within the main part of the profile 



reaches the radius 



I ~ Ro, 



R < T^^Ro, 
R > ifRi,. 



(12) 
(13) 



Such a spreading in density profile may result in a substan- 
tial increase in the width / ^ Rq of the shell when it becomes 
transparent to radiation. Note that condition ( fT3] l coincides 
with the condition in the third line of which corresponds 
to the case of a photon thin outflow. This coincidence may 
give an impression that photon thin asymptotics does not ex- 
ist for hydrodynamically spreading outflows. Nevertheless di- 
rect integration of (|5]l in numerical simulations shows that the 
photon thin asymptotics is valid even for the shell undergoing 
such spreading. 

All asymptoti c solutions for the op tical depth have been 
considered by Mes zaros et al.l (120021) . except the case of 
pure electron-positron outflow. They derived the photo- 
spheric radius considering the expansion in comoving ref- 
erence fram e. Notice, that the photon thin asymptotics is 
obtained in iMeszaros et al.l (|2002|) by assuming hydrody- 
na mic spreading of the outflow found by the same authors 
in (IMeszaros et al.|[T993h . In absence of such spreading (e.g. 
portion of wind considered above) such asymptotics cannot 
be obtained this way. Finally, Tomaet al. (2011) discuss all 
asymptotic solutions, applying them to a relativistic wind. 

We conclude that interpretation of the formula (|7]l it terms 
of photon thick and photon thin conditions given above pro- 
vides additional physical insight to the consideration of opti- 
cal depth in shell and wind models. 

4. PHOTON DIFFUSION 

Up to now we assumed that photons escape the expanding 
outflow when its optical depth decreases to unity. We have 
distinguished two possibilities: 

• in the photon thick case electron number density de- 
creases along the photon path so rapidly that the 
medium becomes too rarified to sustain collisions. 
Most photons however still remain inside the outflow 
after decoupling; 

• in the photon thin case the variation of the electron 
number density along the photon path can be neglected, 
but the mean free path of photons increases with expan- 
sion and eventually exceeds the radial thickness of the 
outflow. 

In all the previous discussion we explicitly neglected the ef- 
fect of radiative diffusion out of the outflow by random walks 
of photons. Diffusion time is given by c = l^/Dc, where 
Ic = Tl is the comoving radial thickness of the outflow, and 
diffusion coefficient is Dc = (cAc)/3 = c/(3crnc), where Ac and 
lie are comoving mean free path of photons and comoving 
electron number density, respectively. 

In order to determine at which radii diffusion becomes im- 
portant one has to compare this diffusion time with comoving 
expansion time of the outflow tc = R/(cT). Taking into ac- 
count (O and dSj we obtain that it happens when the outflow 



1/3 



7.2 X 10'^(£'54Z8?72)'''^ cm, 
5.0x 10"(i'5oAff?72)'''^ cm. 



Rd = (toV^RoI^) 

(14) 

This diffusion radius turns out to be always larger than the 
photospheric radius of photon thick outflows, Rd S> Rph, so 
that diffusion is irrelevant for their description. In the oppo- 
site case of photon thin outflows instead the diffusion radius 
is always smaller than the photospheric radius Rq <C Rpi,- In 
this case most radiation leaves the photon thin outflow not at 
its photospheric radius, but before it reaches the diffusion ra- 
dius, when the outflow is still opaque. In other words, the de- 
coupling of photons from the outflow occurs not locally, as in 
the photon thick case, but near its boundaries where photons 
are transferred to by diffusion. In this sense the characteristic 
radius of the photospheric emission is not the photospheric ra- 
dius found from (jTjl, but the radius of diffusion (fT4l i. Besides, 
the comoving temperature of escaping radiation is different 
from that discussed in Section |3] 

In what follows we consider decoupling of photons from 
photon thick and photon thin outflows separately. 

5. PHOTOSPHERIC EMISSION FROM PHOTON THICK 
OUTFLOWS 

5.1. Geometry and dynamics of the photosphere 

Unlike traditional static sources usually dealt with in astro- 
physics, relativistic outflows may have strongly time-varying 
photospheres. For the portion of wind the optical depth can 
be calculated analytically both at acceleration and coasting 
phases for photon thin and photon thick outflows. The result 



Tir,0,t) = ToRo 



•■sin 6 



-tan 



cr + rcos^ 



0m 



r v/(cr + rcos6')2 + (rsin6')2 



[(cr + rcos6l)2 + (rsin6')2] 



-13/2 



(15) 



where T is the time interval during which photon remains in- 
side the outflow, determined by the equations of motion of 
the photon and of the outflow, and f5,„ = 1-1 /(2Tf). For a 
given laboratory time t the photosphere geometry r = r{9) is 
obtained by equating ( fTSl ) to unity. Then formula ( fTSl l gives 
complete information on the dynamics and geometry of the 
photosphere of portion of ultrarelativistic wind. In order to 
understand this dynamics it is instructive to consider its limit- 
ing cases. 

Firstly, the photosphere of the coasting infinitely lo ng rel- 
ativist ic wind with F = const analyzed by Abramowicz et alj 
([1991 may be recovered from ( fTSl ) with T ^ oo. In that case 
the last term in (fTSl l can be neglected and we have (see e.g. 
IPe^ (12001 -) 



Ro 



= tq 



smf 



-A, 



(16) 



which is a static surface having concave shape, see Fig. |6] 

Secondly, the photosphere of the accelerating infinite wind 
may be obtained from ( fTSl ) for T — > oo and — > oo. It results 
in a qubic equation describing a static surface with curvature 
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Fig. 6. — The shape of photospheres of infinitely long coasting (blue solid 
curve) and accelerating (black dotted curve) winds for Fp/, = 100. Dashed 
line shows the relativistic beaming angle. 




larger than that of the coasting wind, see Fig. |6] In both cases 
these photospheres appear for a distant observer as static spots 
with radius 

P=TTToRo (17) 

and brightness decreasing from the center to the edge. 

Now consider dynamic properties of the photosphere of 
photon thick outflow described by ( fTsT i as seen by a dis- 
tant observer. The arrival time of radiation is defined as 
ta = t-rcosd/c. The equitemporal surface (EQTS) of the pho- 
tospheric emission represents a part of the phot osphere visi- 
ble at a given instant of arrival time ?„, see e.g. iRuffini et alj 
(|2001). We will refer to that surface as Photospheric EQTS 
(PhE). 

PhE of the photon thick outflow has concave shape, see Fig. 
|2]and Fig. [8]for accelerating and coasting cases, respectively. 
This concave PhE in both cases approaches the photosphere 
of infinitely long wind. In the coasting case the approach to 
that surface is only asymptotic, while in the accelerating case 
the photosphere actually reaches it at finite arrival time. The 
external boundary of the PhE for a given ta shown in Fig. 
[8] is defined by the condition that the optical depth for pho- 
tons emitted from the outermost layer of the outflow equals 
unity. Notice that this boundary is wider than the relativistic 
beaming surface (these are tube and cone for accelerating and 
coasting outflows, respectively). As soon as the innermost 
part of the outflow reaches the photospheric radius, i.e. ob- 
server sees the switching off of the wind, the inner boundary 
of the PhE expands with tg. The surface of these boundaries 
is given by ( fT6b in the case of coasting photon thick outflow. 

5.2. Observed flux and spectrum 

The basis of spectrum and flux calculation is radiative trans- 
fer equation for specific in tensity along the ray (see e.g. 
iRvbicki & Lightnianl(IT979h . p. 11) 



dl^_ . _ 

as 



(18) 



where is monochromatic emission coefficient, k^, is ab- 
sorption coefficient and s is distance, measured along the ray. 

Spectral intensity of radiation at infinity on a ray coming to 
observer at some arrival time ta is given by formal solution of 



Fig. 7. — PhE of the photon thick accelerating outflow at different arrival 
times, and the probability density function of photon last scattering integrated 
over angles P(r) (red thick curve). Thick black curve represents the photo- 
sphere of infinitely long accelerating wind. PhEs are illustrated for several 
arrival times with logarithmic spacing by thin black curves. The surface 
cos 9 = (9 is given by p = So and it is shown by thick black line. Dashed 
curves illustrate the maximal visible size p for several arrival times with log- 
arithmic spacing. The PhE at that arrival times is a part of the wind photo- 
sphere limited by the corresponding curves. Range of observed temperature 
of emission under the asymptotic photosphere is illustrated by color, see leg- 
end. HereRp/, = 100i?o- 




0.6 



Fig. 8. — The same as in Fig. |2]for a photon thick coasting outflow. PhEs 
are illustrated by thin curves for several arrival times with linear spacing. 
Thick black curve bounding PhEs correspond to the position of maximal vis- 
ible angles at given ta . Notice that these angles exceed the relativistic beam- 
ing angle, shown by dashed black line. Righ t ed ge of the colored area is the 
photosphere of infinitely long coasting wind (16). Here F = 100. 

this equation 

tuii^, P,ta) = j lAv, r, 9, f )^ {exp[-r(z/, r, 0, f )]} ds 

IAiy,r,9,t)exp[-T(i^,r,9,t)]dT, (19) 
where I^(r,9,t) is the source function, equal to the ratio 
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of emission and absorption coefficients Xy = ju/ Hw, optical 
depth r is an integral of from the point on the ray un- 
der consideration to infinity r = J K^ds, and variables {r,9,t) 
are connected by = t-rcos9 and rsin^ = p, see Fig. [T] 
We use Thomson scattering cross section in comoving frame 
Kc = const. 

Total observed flux is an integral over all rays 



(20) 



where Afi is the solid angle of the observer's detector as seen 
from the outflow in the laboratory frame and Inpdp is an ele- 
ment of area in the plane of the sky. 

In what follows we assume that emissivity is thermal and 
isotropic in comoving frame of the outflow. The laboratory 
source function is then 



2h 



exp 



r(l-ffcos6 



(21) 



■1 



This approximation is justified when the radiation field is 
tightly coupled to the matter. The photospheric emission 
comes from entire volume of the outflow, and the computa- 
tio nal method sketche d above is closely related to that used 
in (iBeloborodovl 120 1 11) where the concept of "fuzzy photo- 
sphere" was introduced. This method will be referred to as 
fuzzy photosphere approximation. 

Most of energy reaching observer is emitted from the region 
near the PhE, where the probability density function along the 
ray 

P(r,0,t) = Po^sxp[-T(r,e,t)] (22) 
ds 

with Po being normalization, reaches the maximum. For this 
reason the dynamics of PhE studied in the previous section 
determines both light curves and spectra of observed photo- 
spheric emission. When the time dependence in this equation 
is discarded this P(r, 6) coincides with the probab ility density 
function of the last scattering defined in (|Pe'erll2008l) . As- 
suming that all the energy comes from the PhE only, i.e. a 
surface instead of the volume discussed above, the computa- 
tion may be reduced to one dimensional integration by substi- 
tution of the function P with a Dirac delta. Such more crude 
approximation, in contrast to the fuzzy photosphere one, will 
be referred to as sharp photosphere approximation. 

For photon thick outflow the optical depth becomes func- 
tion of r and 6 only and the comoving temperature also de- 
pends only on radius. In this respect the photon thick case 
is similar to the infinite wind. Then the integrand in ( |20] | 
does not depend on time and only limits of integration pro- 
vide time dependence due to motion of the outflow bound- 
aries. The probability density function ( |22] | integrated over 
angles is shown in Figs. |2]and|8]for accelerating and coasting 
photon thick outflows. 

The observed flux of photospheric emission from acceler- 
ating outflow is illustrated in Fig. |9]by thick red curve (fuzzy 
photosphere) and by dotted blue curve (sharp photosphere). 
The characteristic raising and decaying time is in both cases 



5t: 



(Rphc) 



(23) 



There is no simple analytic expression describing full light 
curve, however its decreasing part is close to a power law with 
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Fig. 9. — The light curve of photospheric emission from the photon thick 
accelerating outflow in fuzzy (red curve) and sharp photosphere approxima- 
tions (dashed blue curve). Here Rpi, = IOOSq. and At = 2St, see )23t below. 

index -4.7 and -6.5 within fuzzy and sharp photosphere ap- 
proximations, respectively. As minimal duration of the pho- 
ton thick outflow Af is of order Rq/c, then Af ^ 5t and the 
light curve has almost rectangular shape. 

Such accelerating outflow appears to a distant observer as 
a spot with size p = (7;„-(fa/c)^)'/^, for -Ri)/c <ta<0. As 
soon as the PhE reaches the corresponding accelerating in- 
finitely long wind photosphere at fa = the spot size starts to 
increase almost linearly with time p ~ R^ + cta. Finally, as the 
innermost part of the outflow reaches the photospheric radius 
the spot transforms to a ring with rapidly decreasing width 
and brightness. 

The observed photospheric emission of the coasting photon 
thick outflow results in the flux changing as 



■[l-(tp/taf 



with 



R 



ph 



2rf 



(24) 



(25) 



i.e. increase up to the saturation value F„jax oc L, see the rais- 
ing part of the light curve in Fig. [TOl both in sharp and fuzzy 
photosphere approximations. Radius of the visible spot then 
reaches its maximal size ( [Tt] ). As arrival time exceeds tp + lS.t 
the innermost part of the outflow approaches the wind photo- 
sphere ( fTSI ) along the line of sight and the spot transforms to 
a ring, the flux decreases rapidly in both approximations 



F <^ti 



1 



1 



(fa-Af)2 tl 



(26) 



For ta ^ Af it behaves as F cx f~^, see the decreasing part 
of the light curve in Fig. [TO] Similarly to the accelerating 
photon thick outflow the light curve for Af ^ tp has almost 
rectangular shape due to the fact that its increase and decay 
times are much shorter than Af . 

Accelerating photon thick outflows exhibit photospheric 
spectra close to thermal ones, see Fig. [TT] In ultrarelativis- 
tic case spectra computed using both sharp and fuzzy photo- 
sphere approximations are very close to each other Both have 
small deviations from thermal spectrum in the low energy part 
with the corresponding Band low energy indices a = 0.82 and 
a = 0.71, respectively. 

In contrast, the spectrum of photospheric emission of the 
coasting photon thick outflow is significantly wider that the 
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Fig . 11 . — Instantaneous spectrum of photospheric emission of accelerating 
photon thick outflow in continuous (blue thick curve) and sharp photosphere 
approximation (black thick curve). Dashed red curve represents the thermal 
spectrum with the temperature at the line of sight T^os- Lorentz factor at 
photospheric radius is F,,/, = 100. 

thermal spectrum, see Fig. [121 Low energy part is described 
by a power law with Band indices respectively a = 0.34 and 
a = 0, for sharp and fuzzy photosphere approximations. 

After initial phase of evolution, namely rising of the low- 
energy part, spectra do not evolve until observer detects emis- 
sion from the innermost part of the outflow. At that moment 
there is a transition to another phase characterized by the fast 
decrease of both temperature and flux. Considering time- 
integrated spectrum we find that as characteristic times of the 
first and third phases are much less than that of the second 
one, the spectrum is close to the instantaneous one described 
above. 

5.3. Comparison with previous works 

No w we compare the r esults with obtained by other meth- 
ods. Beloborodovl (1201 ll) considered the photospheric emis- 
sion from infinitely long wind both at acceleration and coast- 
ing phases and solved the corresponding steady radiative 
transfer equation. His main conclusion is that in addition to 
usual relativistic beaming leading to anisotropy of radiation 
in laboratory frame, in the coasting wind another anisotropy 
in the comoving frame of the outflow is developing. This co- 
moving anisotropy results from the fraction of photons which 
already underwent their last scattering in the bulk photon field 



of the outflow. The anisotropy of such photons grows with in- 
creasing radius for geometrical reasons. Since the amount of 
such photons increase with radius the entire photon field be- 
comes increasingly anisotropic. 

For the finite photon thick outflow the radiative transfer 
problem becomes explicitly time dependent. The expanding 
outflow at a given laboratory time spans only a finite part of 
the probability density distribution shown in Figs|7]and|8] that 
results in difference in observed spectra for finite and infinite 
cases. Interesting consequence of ultrarelativistic motion of 
the outflow is that even geometrically thin outflow I R at 
a given arrival time spans large interval of laboratory radius 
Ar = 2rH. 

The effect of ad ditional comoving an isotropy on the source 
function found bv lBeloborodovl ( |201 Ih is actually small. The 
difference between the probability of last scattering (l22l l and 
the distribution of last scattering in a steady wind does not 
exceed several percent. 

Our method is also similar to the one used bv |Pe^(l2008l) 
and bv lPe'er & Rvdd (1201 ll) who described the late-time pho- 
tospheric emission of switching off relativistic wind consider- 
ing the probability density function (l22l i for the last scattering 
of photons. An additional approximation adopted by these 
authors is the possibi l ity to split radial and angular dependen- 
cies. Actually iPe'eii ( 120081) computes not the traditional en- 
ergy flux understood as energy crossing unit area in unit time, 
but photon flux as number of photons crossing unit area in 
unit time. For this reason his decay law for photon flux at late 
times is F"''(ta) cx t~^. Lorentz transformation of the photon 
energy from the comoving frame to the laboratory one results 
in additional multiplier (1 - in the energy flux that leads 
to the observed fl ux F (x t~^, wh ich agrees with our result in 
Eq. (|26ll, see also lPe'er & Ryde (2011). 

We conclude that the fuzzy ph otosphere approximati on in 
fact follows closely methods of iPe'er & Rvdd (1201 Ih and 
iBeloborodovl (1201 11) . In fact we obtained similar results for 
the probability of last scatt ering as more sophis ticated treat- 
ment of radiation transfer (IBeloborodovl 120 i T) . The sharp 
photosphere approximation provides good description of light 
curves including their raising and decaying parts. The ob- 
served spectrum from accelerating outflow is also well de- 
scribed in this approximation, while there is some difference 
for the coasting case. The advantage of sharp photosphere ap- 
proximation for computing observed light curves and spectra 
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Fig. 13. — Evolution of PhE for the photon thin coasting outflow and di- 
mensionless radiative diffusion flux, corrected for adiabatic cooling (red thick 
curve). PhEs shown by thin curves correspond from left to right to anival 
times ta = (ff /5,2ff /5, . . . ), see j29t . Thick black curve bounding PhEs cor- 
respond to the surface cos 9 = /3. Relevant range of observed temperature of 
photospheric emission is illustrated by color, see legend. Here F = 100. 



is evident for intrinsically variable and dynamic outflows. 

6. PHOTOSPHERIC EMISSION FROM PHOTON THIN 
OUTFLOWS 

Now we turn to photon thin outflows. In Section |4] we 
pointed out that most of the radiation leaves the outflow not 
at its photospheric radius, but earlier, before the diffusion ra- 
dius. Given that opacity of the outflow is still large, the emis- 
sion escapes only from a very narrow region near the outer 
boundary of the outflow. The probability density function ( |22] | 
is strongly peaked there and the photospheric emission for a 
given arrival time originates from this narrow region. Sharp 
photosphere approximation is thus completely justified in this 
case. 

Since all the radiation is emitted from the PhE we briefly 
discuss its geometry and dynamics. The PhE of the photon 
thin outflow is similar to EQTS of infinitesimally thin con- 
stantly emitting rel ativistic shell cons idered firstly bv' Couder^ 
(11939) and then bv lReesI ([19661 [19671) . The EQTS of this shell 
appears to a distant observer as an ellipsoid with axes ratio 
equal to F. However the PhE of photon thin outflow is not the 
entire ellipsoid: it is only a part of that surface, see Fig. [T3] 
The external boundary of the PhE for a given is defined by 
the condition that photons emitted from the outermost layer 
of the outflow toward observer leaves the outflow. In the pho- 
ton thin asymptotics this surface coincides with the relativistic 
beaming cone. 

We again start with the radiative transfer equation ([TSl l. In 
contrast with the photon thick case, here the source function I 
in ( [T9] l strongly depends on both r and t. The main process by 
which photons are coupled to the matter is Compton scatter- 
ing which conserves the number of photons. Since opacity is 
large other processes which do not conserve the photon num- 
ber lead to local thermodynamic equilibrium with thermal co- 
moving radiation intensity 1^, number density and spectrum 
of photons in the outflow. Hence we us e the Rosseland radia- 
tive diffusion approximation (see e.g. i Rvbicki & Lightnianl 
([T979.) . pp. 39-42), that we derive from the radiative transfer 



equations ( [TsT i for expanding outflows. Details of calculations 
are presented in Appendix. 

It is useful to introduce the function Lc(^, t) = (t / to)^^^Ic{Ci 
which accounts for the adiabatic cooling of radiation in ex- 
panding outflow. Here both ^ and time t are measured in lab- 
oratory frame, while 4(^,0 is measured in comoving frame. 
The diffusion equation wi th these variables is w ell behaved in 
ultrarelativistic limit, see iBeloborodovl (1201 ll) . We then ob- 
tain 



dct 3Ro dS.^ 



■0, 



A: 



Notice that the diffusion coefficient is explicitly time depen- 
dent due to the expansion of the outflow. In Appendix we also 
discuss initial and boundary conditions and obtain an approxi- 
mate analytic solution for the radiation field inside the photon 
thin outflow, which is represented in Fig. [13] The raising part 
of the corresponding flux of Lc through the external bound- 
ary of the outflow scales as f'/^, while its decaying part is 
quasi-exponential one. Consequently, while the diffusion in a 
static object gives the flux decreases as f"'''^, in our case the 
observed flux (l20i ) is a more slowly decreasing function 



Fact 



-1/6 



up to arrival time of diffusion 



2rf 



(28) 



(29) 



where large part of energy has left the outflow already. At 
this moment the energy decrease due to diffusion becomes 
substantial even in the deepest parts of the outflow and later 
the observed flux decreases quasi-exponentially with arrival 
time. 

The comoving temperature of radiation on the photosphere 
is determined by the balance between the energy diffusion 
from the interior of the outflow and radiative losses and it is 
much smaller than the temperature in the interior. The vari- 
ation of observed temperature across the PhE is small, see 
Fig. [T3]and hence the observed instantaneous spectrum is very 
close to the thermal one and peaks near the observed temper- 
ature on the line of sight. We find that the latter decreases as 

ta ^^^^^^ in contrast with adiabatic law ta^^^. However, at diffu- 
sion radius both temperatures coincide giving for the line of 
sight temperature 



Tlos - keV. 



(30) 



The time integrated spectrum has a Band shape with a cut-off 
near the temperature of transition from acceleration to coast- 
ing, see Fig. [T4| Low energy part of the spectrum has the 
slope a=\, while high-energy part has /? ~ -3.5. 

7. DISCUSSION 



The expression ( |29l ) gives an estimate for duration of pho- 
tospheric emission of photon thin outflows. When available 
observed spectra are integrated on time intervals comparable 
to ( |29] | the observed spectrum of photospheric emission is ex- 
pected to have Band shape. Thus, starting from comoving 
thermal spectrum for the photospheric emission we obtain for 
the first time an observed spectrum which may be well de- 
scribed by the Band function with high energy power law in- 
dex (3 being determined by the density profile of the outflow. 
We find this result quite remarkable. 
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Fig. 14. — Time-integrated spectrum of pliotosplieric emission of photon 
thin outflow (thick curve, rj = 100, Rd = 10^i?o), superimposed with two in- 
stantaneous spectra of that emission, corresponding to arrival time of photons 
emitted at the moment of transition from acceleration to coasting (dashed 
green curve on the right) and to arrival time of diffusion (dashed red curve on 
the left). 

Notice that non thermal spectra as a result of convolution of 
thermal ones ov er time ha s been introduced for afterglows of 
GRBs by Blinn ikov et al.l (Il999h . Double convolution over 
EQTS and arrival time is also one of the key ideas in the 
fireshell model (Ruffini et al. 2003). 

Band spectra in photospheric models of GRBs have been 
obtained by now only assuming additio nal dis sipation mecha- 
nisms such as m agnetic reconnectiori (iGiann ios 2006), col- 
lisional he ating ("Beloborodov '2010") and internal shocks 
jToma et a l. 201 1; Rvde et al. 201 lb). In our model such ad- 
ditional assumption is not required. 

It is even more remarkable that GRBs appear to be the only 
known objects in nature able to reach the photon thin asymp- 
totics in their ultrarelativistic expansion. For thermally ac- 
celerated relativistic plasmas which are discussed in connec- 
tion with their po ssible synthesis in ground based laboratories 
(e.g. iRuffini etal . (2010)) it is unreachable. The photon thin 
asymptotics is reached if the optical depth (O 

ro»4r4±=4xlO^^^. (31) 

GRBs clearly can satisfy this constraint as the contribution of 
baryons tq ~ 3.5 x 10'^£'54r72'^8''8'- 

We obtained both time integrated and instantaneous ob- 
served spectra of accelerating photon thick outflow which are 
close to thermal one with sm all deviations in the Rayleigh - 
Jeans part, in agreement with Grimsrud & Wasserma"nl(ll998l) . 

Time integrated observed spectrum of coasting photon thick 
outflow is broader than thermal one, and deviates from it both 
in low- and high -energy parts. This broadening is also found 
by iBeloborodovl (|2010) using Monte-Carlo simulations, but 
our spectral index of the low energy part is a = instead of 
his a ~ 0.4. This difference may be attributed to the fact that 
Monte-Carlo simulations do not account for stimulated emis- 
sion. Interestingly, we found a ~ 0.34 in the sharp photo- 
sphere approximation. 

As discussed earlier in Section [3] each photon thick out- 
flow always contains a photon thin layer with depth ^,/„„ = 
(R + 1) /(2rf) located near its outer boundary. Radiation dif- 
fused out from this part of the outflow arrives to observer 
first and modifies the initial part of the light curve and the 
corresponding spectrum of the outflow. The diffusion length 



1 /2 

£.D = [(R + l)^ /{v^toRq)\ remains always within this pho- 
ton thin layer ^,/„„ and our solution for photon thin outflow 
is applicable for description of this early emission. The pho- 
tospheric emission from well developed photon thick outflow 
and its late time behavior occurring when the outflow crossed 
the photospheric radius may be described either by fuzzy or 
by sharp photosphere approximations. 

When the outflow becomes transparent in the transition 
from photon thick to photon thin conditions, the observed 
time integrated spectrum will contain both Band component 
produced by the early emission from the photon thin layer and 
thermal-like component coming from the photon thick part 
superimp osed. This may be the reason why in most GRBs an- 
alyzed bv lRvde & Pe'erl (|2009) there are both power law and 
black body components. 

We presented analytic expressions for the photon flux in 
the simple model of the portion of wind. With more complex 
density profile composed of presumably many shells the light 
curve is expected to be variable and arbitrarily complex. The 
minimal time scale of variability is given Eq. (|29] l and it may 
be very small for small baryon loading. It is necessary to 
emphasize that the decaying part of the light curve follows 
for photon thick outflows. Steeper decay of the light curve of 
photospheric emission it is a clear signature of the photon thin 
outflow. 

The photospheric emission should be additionally identified 
by the spectral analysis. In particular, power law spectra ex- 
tending to high energies above 10 MeV cannot be produced 
by the photospheric emission unless additional mechanisms 
are involved. What we have shown here, though, is that the 
observed spectrum may not necessarily be close to the thermal 
one. 

8. CONCLUSIONS 

In summary, in this paper we proposed a unified treatment 
of ultrarelativistic outflows, which originate both from instan- 
taneous and from continuous energy release, with respect to 
photospheric emission. We have to stress that these two cases 
are discussed in the literature and are respectively referred to 
as shells and winds. Instead of this traditional division we 
propose a new physically motivated classification, which in 
our opinion helps to understand in particular why geometri- 
cally thin shell may appear as thick wind with respect to the 
photospheric emission. For this reason we re-examined the 
existing scattered literature and pointed out the advantage of 
the proposed classification. 

We studied geometry of photospheres in generic relativistic 
outflows. As we are interested in appearance of the photo- 
sphere to a distant observer, we introduced the notion of pho- 
tospheric equitemporal surface and described its dynamics. 

We computed both energy flux and observed spectra of pho- 
ton thick outflows in two approximations, derived from the 
radiative transfer equation. In our fuzzy photosphere approx- 
imation the effect of simultaneous emission from entire vol- 
ume of the outflow is taken into account. We also used another 
computationally more simple sharp photosphere approxima- 
tion which is shown to reproduces well both light curves and 
spectra. These results generalize the corresponding results in 
the literature for steady relativistic winds. 

In photon thin outflows most of radiation is shown to orig- 
inate not at its photospheric radius, but at smaller radii due 
to radiation diffusion. Starting from the radiative transfer 
equation for time dependent outflows we derived the diffu- 
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sion equation and obtained approximate analytic solution for 
the energy flux. We present both instantaneous and time inte- 
grated observed spectra. The latter are well described by the 
Band function. For our simple density profile we find values 
for the low energy power law index a=l and the high-energy 
power law index /? ~ -3.5. 
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APPENDIX 

For large opacities the distribution function of photo ns in comoving reference frame is close to isotropic one and the radiative 
diffusion approximation is accurate. Following Belo borodovl (1201 ll) we use spectral intensity in comoving frame J^(t,^,ij). 
Starting from the radiative transfer equation ( fTsl ) along the ray s in laboratory frame, we transform all variables except time t, 
depth ^, and distance s into comoving reference frame 



Theory of photospheric emission from relativistic outflows 

and integrating over comoving frequency v we have 

IdJ fj. dJ l-i/ dJ ,r/3 1-/1- 



c dt TVdS, vt-idn 



+ 4- 



V vt-S, V 
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(33) 



where fi = cosOc, 9c is the photon angle with respect to the radial direction in comoving frame, V = r(l + is Doppler factor, 
J = f Judv is total photon intensity, 5 = J S,jdv is total source function, k = J Ki,Ji,dv is effective opacity, n^, = an^ is opacity 
in comoving frame. 
In the case of small deviations from isotropy decomposition 



7 = 7o(f,0 + M-/i(f,0 



(34) 



could be applied. Introducing it into ( |33] | and integrating it over Vdn and over Vjidji after some algebra for coherent scattering 
with 5 = 5o = /o we have 
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Diffusion approximation is based on slow variation of total flux through the entire sphere Li = Ji(t //o)^ over mean free path, so 
that ^ = 0, and it provide Ji from the equation ( |36] |. Inserting this into (l35T l after simple but tedious calculations in ultrarelativistic 
/3 ~ 1 photon thin case F^^ ^ vt for function L = Joit /U)fl^ we obtain the diffusion equation 



da~ 3Ro dS} ' 



^0, 



A: 



(37) 



This equation should be supplemented with boundary conditi ons. There are two types of boundary conditions used frequently: 
free-streaming, for example in two-stream approximation (Ryb icki & LightmanI (Il979h . pp. 42-45), and ze ro boundary condi- 
tions, that can be used as replacement for free-streaming for "extrapolated boundary" (iHaskell et al.lll994 ). We find that the 

position of "extrapolated boundary" ^ = -k^-j^ (k is a constant of order unity, dependent on the approximation used for free- 
streaming description) for the main part of emission is very close to the real boundary, and in the case of zero boundary conditions 
L|^=o = L\^=i = there is a series expansion of solution, that for initial conditions L(^,fo) = 1 gives 
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exp 
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(38) 



This solution in comparison with numerical one with free-streaming boundary conditions is accurate to a few percents. 
The flux of L is characterized by an initial burst and then tends to the asymptotic solution, that corresponds to fo = 0, with flux 



3.2 



F(t): 



■d2 



0,exp 



4A7r2c^\3 



where ^2 is the Jacobi elUptic theta function, see Fig. [T3] The peak of the flux of L is near the diffusion time 



to- 
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1/3 



(39) 



(40) 



and "extrapolated boundary" ^ = -klilA/Ro)^'^ <C I is very close to real one as A <C 1, that ensures the accuracy of (|38T l. 



